Hole subband mixing and polarization of luminescence from quantum dashes: a 

simple model 
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In this paper, we address the problem of luminescence polarization in the case of nanostructures 
characterized by an in-plane shape asymmetry. We develop a simple semi-qualitative model revealing 
the mechanism that accounts for the selective polarization properties of such structures. It shows 
that they are not a straightforward consequence of the geometry but are related to it via valence 
subband mixing. Our model allows us to predict the degree of polarization (DOP) dependence on 
the in-plane dimensions of investigated structures assuming a predominantly heavy hole character 
of the valence band states, simplifying the shape of confining potential and neglecting the influence 
of the out-of plane dimension. The energy dependence modeling reveals the importance of different 
excited states in subsequent spectral ranges leading to non-monotonic character of the DOP. The 
modeling results show good agreement with the experimental data for an ensemble of InAs/InP 
quantum dashes for a set of realistc parameters with the heavy-light hole states separation being 
the only adjustable one. All characteristic features are reproduced in the framework of the proposed 
model and their origin can be well explained and understood. We also make some further predictions 
about the influence of both the internal characteristics of the nanostructures (e.g. height) and the 
external conditions (excitation power, temperature) on the overall DOP. 

PACS numbers: 78.67.Hc,73.22.-f 



Quantum dashes (QDashes) [iHal are quai-zero- 
dimensional nanostructures characterized by strong in- 
plane asymmetry. TypicaUy, their width is on the or- 
der of several to a few tens of nanometers, while their 
length may be on the order of hundreds of nanometers 
[l|-[j]. The electronic structure and optical properties of 
these structures have been investigated both experimen- 
tally [3-[al and theoretically [3, [ol-Ulj , which is motivated 
by their favorable properties from the point of view of 
photonic applications at telecom wave lengths Q. For 
those and other applications, the understanding of po- 
larization properties (expected to be anisotropic due to 
shape asymmetry) is crucial. 

In Fig. (T^a) we show the polarization-dependent pho- 
toluminescence spectrum obtained at 300 K from an 
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FIG. 1: (a) Room temperature photoluminescence spectra 
of an ensemble of QDashes at two orthogonal linear polar- 
izations (red solid line: parallel to the structure elongation, 
blue dashed line: perpendicular to the structure elongation), 
(b) The corresponding degree of polarization. The apparent 
structure in the high-energy part results from irrelevant fluc- 
tuations due to water vapor absorption in the aera of weak 
photoluminescence signal. 



ensemble of epitaxially grown self-assembled InAs/InP 
QDashes preferentially elongated in the [110] direction, 
similar to those studied in Ref. [2, S]. Fig. [TJb) shows 
the corresponding degree of polarization (DOP), 
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where ^| | , ^_l are the intensities of luminescence compo- 
nents polarized along and perpendicular to the QDash 
elongation axis. The anisotropic shape of the QDashes 
leads to a high DOP of the emitted radiation (see also 
Ref. [8|). This polarization cannot be a trivial conse- 
quence of the confinement shape at the level of a single- 
band effective mass and envelope function description as 
the emission is governed by the interband dipole moment 
that depends on the Bloch parts of the wave functions 
and not on their envelopes. Therefore, a more general 
description is needed in order to take subband mixing 
into account. Quantitatively exact modeling of such non- 
trivial optical properties of nanostructures must be based 
on atomistic or multi-band kp methods [10, II4I . How- 
ever, it may be useful and interesting to have a simple 
semi-quantiative model that would not only elucidate the 
physical mechanism of the polarization but also yield an 
estimate of the expected effect in terms of the shape pa- 
rameters. 

In this contribution, we present a "minimal" theory 
that is able to account for the observed DOP in the 
QDash luminescence. First, we show that the polariza- 
tion is proportional to the degree of heavy- light hole mix- 
ing. Then we estimate the degree of the latter within a 
semi-quantitative model and show that it is related to 
the shape asymmetry of the structure. We provide a 
very simple formula relating the DOP to the QDash di- 
mensions. 



The polarization of light emitted in a recombination 
process from a conduction band (cb) state ^c(^) to a 
valence band (vb) state ^v('") is determined by the in- 
terband matrix element of the dipole moment operator 
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Here we have performed the standard separation of 
length scales in the integration, dx is the bulk interband 
dipole moment between the initial cb (say, with spin up) 
and the final valence subband, and 

ax= cPril}l{r)'4>^{r) 

is the envelope function overlap, where we denoted the cb 
envelope wave function by ipcif) and the components of 
the vb envelope by ipxlr), where A is the subband index. 
For the initial state in the spin-up conduction subband, 
only two transitions yield non- vanishing in-plane compo- 
nents of the dipole moment: to the -1-3/2 heavy hole band 
and to the —1/2 light hole band. The corresponding bulk 
matrix elements are 
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where the two components of the in-plane vectors refer 
to the Cartesian coordinates x, y bound to the crystallo- 
graphic axes. 

In order to study the properties of a nanostructure 
elongated along the [110] direction we define the unit 
vectors along and transverse to the elongation direction 



ei 



G-x 67, 



V2 ' 



et 



V2 



The components of the interband dipole moment along 
these directions are 
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Note that we use the standard definition of the basis 
functions [13l . Il4| (which differs from that used in many 
papers employing the kp theory p^. [la|). 

The intensities of the linearly polarized components 
of the emitted radiation along and transverse to the 
structure are proportional to |dip and |dtP, respectively. 
Hence, the DOP is P = {\dy\^-\dt\^)l{\dy\^ + \dt\^). From 
Eqs. ^ and (O, one finds 
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From these equations, it is clear that mixing between 
heavy and light hole contributions in the confined hole 
state can lead to anisotropy of emission polarization with 
respect to the structure geometry, depending on the rel- 
ative phase of the light-hole and heavy-hole components. 
A more quantitative conclusion may be achieved if one 
assumes that the lowest hole state is predominantly of 
heavy hole character, with an admixture of light hole 
states. This is justified in many structures since the light 
hole states are shifted in energy with respect to the heavy 
hole states due to confinement and strain and the inter- 
band coupling elements are relatively small. It follows 
from the structure of the Kane hamiltonian [ij, [l^ that 
the heavy hole state with the angular momentum -1-3/2 
is coupled in the leading order of perturbation to both 
light hole states. However, the coupling term between 
this heavy hole state and the -1-1/2 state is proportional 
to kz which means that only states excited along the z 
direction are coupled. Due to the strong confinement in 
the growth direction these states have a very high energy 
and their contribution is expected to be small. Thus, the 
hole state is essentially composed of the -1-3/2 hh com- 
ponent with an admixture of a —1/2 Ih component. In 
the 1st order of perturbation theory one finds for the Ih 
admixture 



V^-l/2(^) 



1 



1 



dV'0:(r')T^V'3/2(r') 



AK 



Ih 



V'-03/2(t-), 



where {(/),i} is any complete set of functions, V is the 
relevant inter-subband element of the kp Hamiltonian, 
A£^ih — E\^ — E\ is the energy separation between the 
heavy and light hole subbands, the differences between 
energies of various Ih states have been neglected, and the 
last equality is obtained from the completeness relation 
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With our choice of basis states, the element of the Kane 
Hamiltonian coupling the relevant states is [ij] 

where kj = id/dxj, mp is the free electron mass and 7^ 
are Luttinger parameters. Only the imaginary part of V 

contributes to Im 
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where we used the components relative to the structure 
elongation. 
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FIG. 2: (a) Theoretical results for the room temperature pho- 
toluminescence spectra of an ensemble of QDashes at two or- 
thogonal linear polarizations (red solid line: parallel to the 
structure elongation, blue dashed line: perpendicular to the 
structure elongation), (b) The corresponding degree of polar- 
ization. The energies are shown relative to the ground state 
at the majcimum of the distribution. 



In order to obtain a general estimate without relying 
on specific information on the confined wave functions we 
note that the dominating contribution to luminescence 
originates from electron and hole wave functions with 
the same quantum numbers and neglect the differences 
between their exact shape, so that 03/2 ~ 1. Assuming 
the simplest rectangular box model for the confinement 
we find 
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where L, D are the dash length and width, respectively 
and ni, nt are the corresponding quantum numbers in the 
rectangular box model. Hence, the DOP is 
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The formula ([3]) is the main theoretical result of our 
study. In spite of its simplicity and approximate char- 
acter, it reproduces experimental data reasonably, as 
can be seen by comparing the measured DOP shown in 
Fig. [T] with the results of theoretical modeling presented 
in Fig. [5] . Here, we plot the results obtained from an 
ensemble of QDashes modeled by a set of boxes (with in- 
finite potential walls) characterized by the fixed height to 
width ratio of HjD = 1/3, identical length L = 150 nm 
and a gaussian distribution of ground state energies with 
the standard deviation of 44 nieV (which uniquely deter- 
mines the size distribution) . This set of parameters cor- 
responds to the geometry of the QDash ensemble [3, 0| 
and to the inhomogeneous broadening of the low temper- 
ature photoluminescence spectra (not shown here). We 
have chosen the effective masses of electrons and holes 
to be TOe = O.OTto-o and TOh = O.Stoq, respectively, where 
Too is the free electron mass, which yields the lowest op- 
tical transition corresponding to a bright state excited 



in the transverse direction (ni = l,nt = 2) at the en- 
ergy IS.E = 41 meV above the ground state. Fitting to 
the overall polarization value observed in the experiment 
yields the heavy-light hole separation ASih — 30 meV, 
which is a reasonable value (actually, however, this en- 
ergy should be correlated with the dot height, which we 
neglect). The results correspond to the thermal distribu- 
tion of electron and hole occupations at 300 K, assuming 
weak enough excitation to keep the Fermi system non- 
degenerate. 

While the overall DOP follows from the value of the 
heavy-light hole energy splitting which is a fitting pa- 
rameter, the agreement of the features appearing in the 
energy dependence is remarkable. Both in the theory and 
experiment, the DOP is roughly constant in the area of 
the maximum luminescence (820 meV in the experimen- 
tal results). Slightly to the higher energies, the DOP 
drops. This can be understood from Eq. ([3]): At higher 
energies, the states excited along the structure {n\ > 1) 
have a larger contribution, which reduces the polariza- 
tion. Farther towards the higher energies, the state ex- 
cited perpendicular to the structure (rit = 2) starts to 
contribute strongly. At 300 K, the occupation of this 
state is non-negligible. According to Eq. ([3]), for low val- 
ues of ni, the DOP of the luminescence from this state 
is four times higher than that of the ground state. This 
leads to a clear maximum in the DOP observed in the 
experiment and reproduced in the model. On the low 
energy side of the luminescence maximum, the theory 
predicts a decrease in the DOP due to the increasing con- 
tribution of larger dots, while the ground state emission 
dominates. The discrepancy with respect to the experi- 
mental value is due to the abrupt decrease of the detec- 
tion sensitivity below 0.78 eV (InGaAs detector) causing 
the increased uncertainty of the experimental results in 
the area of weak luminescence signal. 

In general, Eq. ^ yields many predictions for the DOP 
dependence on various system parameters. For struc- 
tures with a small asymmetry {D ^ L), the DOP should 
increase with growing L. However, this increase satu- 
rates and for strongly elongated structures (L :$> D) at 
low temperatures and weak excitations, where the ground 
state emission dominates, one finds P ^ 1/D'^, inde- 
pendent of L. The DOP should decrease and then in- 
crease with both temperature and excitation power as 
both these factors lead to an increased contribution from 
the excited states, first along the structure (which de- 
creases the DOP), then transverse to the structure (which 
increases the DOP). Finally, semi-speculatively and be- 
yond the presented formal model, one can expect that 
taller structures (larger H) should generally have smaller 
heavy-light hole separations and therefore should show 
stronger polarization. 

In conclusion, we have formulated a semi-quantitative 
model that explains the polarization of light emitted by a 
nanostructure with a strong in-plane asymmetry (a quan- 
tum dash) by relating it to the hole subband mixing. The 
model yields a simple estimate of the DOP in terms of 



the geometrical and material parameters. We have tested 
our model against a room temperature measurement per- 
formed on an ensemble of InAs/InP QDashes. We have 
shown that the model not only accounts for the overall 
DOP for a reasonable value of the assumed heavy-hole 
splitting but also reproduces the features observed in the 
energy dependence of the DOP and allows us to inter- 
pret them in terms of the contribution to the lumines- 
cence from various excited states. Based on our model, 
we have made further predictions for the dependence of 
the DOP on the structure parameters and experimental 



conditions. 
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